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Abstract We consider a belt of small bodies (planetesimals, asteroids, dust parti- 
Q^l cles) around a star, captured in one of the external or 1 : 1 mean-motion resonances 

I ■ with a massive perturber (protoplanet, planet). The objects in the belt collide with 

O I each other. Combining methods of celestial mechanics and statistical physics, we 

> . calculate mean collisional velocities and collisional rates, averaged over the belt. 

^ ' The results are compared to collisional velocities and rates in a similar, but non- 

resonant belt, as predicted by the particle-in-a-box method. It is found that the ef- 
fect of the resonant lock on the velocities is rather small, while on the rates more 



substantial. At low to moderate eccentricities and libration amplitudes of tens of 
degrees, which are typical of many astrophysical applications, the collisional rates 
\^ ' between objects in an external resonance are by about a factor of two higher than 

, those in a similar belt of objects not locked in a resonance. For Trojans under the 

ff^ ' same conditions, the collisional rates may be enhanced by up to an order of mag- 

nitude. The collisional rates increase with the decreasing libration amplitude of 
f-*) ■ the resonant argument, depend on the eccentricity distribution of objects, and vary 

I from one resonance to another. Our results imply, in particular, shorter collisional 

. lifetimes of resonant Kuiper belt objects in the solar system and higher efficiency 

of dust production by resonant planetesimals in debris disks around other stars. 
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1 Introduction 

Resonances between small bodies and giant planets are common in the solar sys- 
tem and, very likely, in planetary systems around other suns. Some of the asteroid 
families in the present-day solar system are locked in various resonances. Exam- 
ples are Greeks and Trojans in the 1:1 resonance with Jupiter, the Hilda group of 
asteroids in a 3:2 resonance with it, or the Koronis family in a 5:2 resonance. In 
the Edgeworth-Kuiper belt, Twotinos and Plutinos reside in the 2:1 and 3:2 reso- 
nances with Neptune, respectively. In the early solar nebula, small planetesimals 
may have been brought to resonance locations b y gas drag and got captured in res- 
onances with protoplanets, too ( Golds ,1975 ; Marzari & Weidenschilling. 2002h . 
Around other stars, resonant belts of planetesimals are thot ight to be re s ponsi- 
ble for clu mpy structure observed in resolved debris disks dWyattl I2003L l2006l : 
iKrivov et a l.. 2007). 

Whether resonant or non-resonant, small bodies in all these systems are sub- 
ject to collisions, whose typical outcome varies from one system to another and 
ranges from perfect agglomeration to full disruption. Accordingly, accurate de- 
termination of collisional velocities and rates amongst the objects is an impor- 
tant ingredient of collisional models. In non-resonant systems, they can be cal- 
culated on the base of the particle-in-a-box approach, in which the system is fol- 
lowed in a reference frame moving around the star with Keplerian circular veloc- 
ity (^pik, 1951; Safronov, 1969; Greenberg et al., 1978, 1991). This method has 
been successfully generalized to include gravitational enhance ment, effects of Ke- 
plerian shear, dy namical friction, and viscous stirring OWetherill & Steward 1 19891 : 
IWetherilll Il990h . It has also been imp lemented as multiannulus schemes to treat 
a ran ge of distances around the star ( Spaute et all Il99ll : IWeidenschilling et"an 
Il997h . Box-based models ha ve been successfully applied for evolution of main- 
belt a steroids fc ampo Bagati n et all 119941: iBottke et alT.ll994 :lDavis & Farinellal 
Il997l : iBottke et al.i. ,2005), Kuiper-belt objects (Sterri [19951 \l99^, debris disks 



of other main-sequence stars (Thebault et al., 20 03), and accumu lation of plan- 
etesimals during planetary system formation (e.g. iLissaued I1993L and references 
therein). 

However, the box method essentially involves an assumption that the angular 
orbital elements of the particles (argument of pericenter and mean anomaly) have 
a uniform distribution. Accordingly, the population of objects in question has a 
rotationally-symmetric structure. This basic assumption is violated in the case of 
a resonance locking. For particles locked in a resonance, the angles vary in such a 
way that a certain combination of them, called resonant argument, librates around 
a certain value with a certain amplitude. As a result, the resonant population may 
be clumpy. One might expect the collisional rates, and possibly relative veloci- 
ties, of particles in such a clumpy structure to systematically deviate from those 
predicted by the particle-in-a box-method. 

The problem can, and has been, studied with direct A'^-body codes. The diffi- 
culty that the "real" number of bodies strongly exc eeds a numerically affordable 
one c an be overcome in the following way (see, e.g.. lThebault & DoressoundiramL 
I2OO3L and references therein). One treats a limited number Nnum of test particles 
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and assigns to these particles an inflated radius such that the total optical depth 
of the numerical sy stem is equal that of the "real" one: NmmiRlum = NreaiR%ai- 
Another possibility dWyatti |2006|) is to use a sort of a "local" particle-in-a-box 
method. An A'^-body code picks up a particle at a location of interest and looks at 
all particles in its vicinity. Then their velocities relative to the target particle and 
their total cross section are calculated. Finally, a particle-in-a box-like approxima- 
tion in that local region gives the collisional rates there. 

The aim of this paper is to make an analytic calculation of collisional rates 
and collisional velocities within a family of objects locked in one of the mean- 
motion resonances with a nearby perturber. Our general method, which com- 
bines classica l celestial mechanics with statistical physics, was published earlier 
dKrivov et al.l [2005, 2006, 2007). The specific approach used here, with a fo- 
cus on non-uniformly distribut ed angles, bears close resemblance to the approach 
by DeirOro and collaborators dDell'Oro & PaolicchiL[l998l : lDell'Oro et alill998L 
l200lLl2002h . 

However, the latter is rather tuned to explore collisions in a finite set of in- 
dividual objects with known orbits (e.g., a given asteroidal family), whereas this 
paper analyzes "infinite" set of fiducial objects with a continuous distribution of 
orbital elements. First, for two arbitrary particles in the resonant family, we for- 
mulate the collisional condition, compute the collisional probability per unit time, 
and calculate the relative velocity of the two particles at the collision point. Sec- 
ond, we derive the distributions of the angular elements of the objects from the 
condition of the resonant libration. Finally, we obtain the average collisional ve- 
locities and rates in the whole ensemble by integrating relevant quantities over all 
possible pairs of colliders, weighted with the "resonant" distributions. 

Section 2 describes the resonant family of objects: essential features of the 
resonant dynamics, simplifying assumptions that we make, choice of variables 
etc. Section 3 deals with a binary collision between two given particles. Here, the 
collisional condition is derived and the velocity of collision is evaluated. Section 
4 defines, in the form of integral formulas, collisional velocities and rates in the 
whole ensemble of objects. These are calculated in Sections 5 and 6, respectively. 
Section 7 addresses the 1:1 resonance case, i.e. Trojans. Possible applications are 
discussed in Section 8. Section 9 contains our conclusions. 

2 Resonant belt 

2.1 System 

We will consider the following system. There are a massive body, which we will 
call planet, orbiting the primary, referred to as a star, and a disk of objects (small 
bodies or dust), orbiting the same primary. Motion of each of the objects or the 
planet in 3D is described by six Keplerian orbital elements 

a , e , i , (i> , Q. , X, (1) 

which stand for the semimajor axis, eccentricity, inclination, argument of pericen- 
ter, longitude of ascending node, and mean longitude, respectively. Instead of A,, 
either the mean anomaly M or the true anomaly 6 can be used. The elements of 
the planet will be marked with a subscript p. 
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2.2 Simplifications 

To keep the problem analytically tractable, we make several major simplifying 
assumptions. 

1. First, as in many other studies, we assume a circular orbit of the planet: 
fip = 0. Apart from reducing the complexity of the problem dramatically, this 
assumption enables comparison with the particle-in-a-box results. For the solar 
system, this assumption is natural. For many other systems where planetesimals 
are expected to be trapped in resonances as a result of planetary migration, it is 
reasonable, too, sin ce diss i pative forces that cause migration tend to circularize 
the planetary orbit ('Wvattl l2003h . If dust particles instead of planetesimals are 
considered, and therefore the mechanism of resonant capture is dust transport by 
dissipative forces rather that planet migr ation, re s onanc e capture is known to be 
most efficient for less eccentric planets. 'O uillenI (|200(?). among others, numeri- 
cally investigated the dependence of the dust capture probability on Cp. She found 
that resonant capture is only possible if (fWp/iVf^,)"'/^^^ ^ 1, where SVfp and £Vf^, 
are the masses of the planet and the star, respectively. Furthermore, even if capture 
occurs, already a low planetary eccentricity of ^ 0.05 smears a clumpy resonant 
structure to a rotationally symmetric ring (Remy Reche, pers. comm). 

2. Further, we confine our analysis to small inclinations: i,ip ^ 10°. One 
reason for that is that orbital inclinations of small bodies in the solar system 
and other planetary systems have this order of magnitude. Another reason is that 
higher inclinations drastically redu ce the probabili t y of resonance capture or make 
the re sonant orbits unstable (e.g., IWisdomL 119831 : bancart et al.L 120031: iGallardol 
l2006h . Furthermore, we will compute collisional velocities by simply assuming 
/ = ip = 0, and take into account corrections due to small non-zero inclinations 
only in the calculations of collisional rates. Thus we can assume a 2-dimensional 
geometry as shown in Fig. [T] and introduce overlined variables, measured with 
respect to the planet: 

00 = (H-Xp, 
A, = A, — Xp. 

3. Next, all particles are assumed to have the same "typical" radius. This 
is justified by the fact that we want to focus on geometrical effects caused by 
resonance locking, rather than on the size distribution effects. 

In fact, this assumption does not imply any loss of generality, as long as we 
consider macroscopic objects, whose dynamics is purely gravitational and there- 
fore independent of their size. It would be, however, a simplifying assumption 
for dust, the motion of which is affected by radiation pressure and is therefore 
size-dependent. 

4. We consider the resonant population only and thus assume all objects to 
have the same, resonant, value of semimajor axis: a = a^es- Thus the semimajor 
axis a is treated as a parameter, not a variable. 

5. Finally, our analysis is confined to external and the primary resonances only: 

^res ^ ^p- 

Additional, less important simplifications, are introduced below. 
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2.3 Mean-motion resonance 

An external mean-motion resonance (MMR) arises when the periods, being the 
inverse of mean motions, of a planet Pp and a particle P are in a rational commen- 
surability: P/Pp = {p + q)/p, where p and q are integers. The MMRs are located 
atUres = Uplip + q) / pf^^ . 

If "particles" are small dust grains rather than planetesimals, Ures is shifted by 
a factor of (1 — P)'/^, where (3 is the radiation pressure to gravity ratio. In this 
paper, we set P to zero. 

For an object residing in the resonance, a certain combination of its orbital 
elements, called the resonant argument, 

^ = {p + q)X - pXp ~ qti) = {p + q)X - qW 

librates around a certain value 4> n, often close to 0°, 180°, o r 270° , with a give n 
amplitude A (Murrav & Dermottl ll999l : iKuchner & Holmanl l2003l : IWvatd l2006h . 
In the case of the primary l:l-resonance (p — 1, q — 0), the resonant argument 
simplifies to 

4) = X 

and librates around the = ±n/3, corresponding to the motion in the vicinity of 
the Lagrange points L4 and L5 . 

Next, the objects captured in a resonance tend to gradually pump their orbital 
eccentricities, up to a certain value e,„ax- Both the libration amplitude A and the 
maximum eccentricity e,„„^ depend, in a complex way, on the mechanism of res- 
onant trapping (Poynting-Robertson effect, planet migration, etc.) Furthermore, 
even for a given trapping mechanism, both quantities depend on a multitude of 
physical parameters: planet mass, order of the resonance, etc. Throughout this 
paper, we assume that both A and e^ax can be "pre-determined" by a dedicated 
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dynamical study of a system of interest, and therefore treat them as, free parame- 
ters. 



2.4 Distributions 

In 2D, four orbital elements fully describe the particle's motion: a, e, 9, CO. As a is 
a parameter, only three remaining elements represent the phase spac e variables. 

W e now introduce the distributions of our variables. Following iKrivov et al.l 
(|2005|) . we denote by i^{x,y^...)dxdy... the fraction of particles in the disk with 
arguments [jc, x + dx] , [}',}' + dy] , ... . The (|)-functions are normalized to unity: 



f '\>{x,y,...)dxdy... = 1. 



The mean longitude X is distributed uniformly, and the distribution of the true 
anomaly 9 follows from the Kepler equation: 



(|)(co,9) = (|)(co,A,) 



39 



: (|)(CO,A,) 



/-2 



a- 



where A, in the right-hand side should be calculated from 9 by means of standard 
formulas of Keplerian motion 

X = (0 + M, (2) 
M = E-esinE, (3) 



E l-e d 

tan- = \ tan-. (4) 

2 V l+e 2 

Within the resonance, where <i>Q — A < <i> < <i>o + A, the distributions of W and 
X (or 9) are not independent. Assuming the distribution of <I> within the libration 
width to be uniform, we obtain 

^^{m,l) = ^^H[^o-A<^<^o+A], (5) 
271 ZA 

where H [cond] is a Heaviside function, which equals one if the evaluation of cond 
returns true and zero if cond rstums false. If needed, H can be replaced by a more 
realistic distribution, e.g. sinusoidal. 

Particles caught in a resonance may have eccentricities between zero and a 
maximum value emo.t- To keep the analysis simple, we assume a uniform distribu- 
tion between the two borders, 

<^e{e;emax) = -^H[e <e„ax]- (6) 

^max 



This simplification, like the one for the resonant argument made above, can always 
be lifted by replacing the Heaviside function with a more realistic distribution. 
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-2-1012 -2-1012 -2-1012 -2-1012 



Fig. 2 Spatial distribution of objects, locked in a 4:3 resonance. The grey scale is arbitrary; the 
darker the grey, the higher the "density" of objects. The star at (0,0) and the planet at (1,0) are 
shown with large and small circles. From left to right: dependence on eccentricity, e = 0. 1, 0.25, 
0.5, and 0.75. From top to bottom: dependence on the libration amplitude, A = O.Iti (18°), 0.37t 
(54°), and Jt (180°). 

Figure |2] gives a visual impression of how the distribution ^ describes a res- 
onance population. Plotted is the distribution of cartesian coordinates ^xy{x,y) 
calculated from (|)g)(co, A.) by the virtue of 

(ifxy{x,y)\dxdy\ = (i?^{m,l)\dmdX\. (7) 



The upper plots, drawn for a small libration width A, show "loopy", pretzel- 
like stru ctures, well-known to be t ypical of the synodic motio n of resonant parti- 
cles (cf. lMurrav & Dermottl[l999l:lKuchner & Holmanll2003h . The middle panels 
show how a larger libration amplitude dithers the distribution. Yet more fuzziness 
will obviously come after convolving Eq. (|5]l with a distribution of eccentricities 
Finally, the lowest panels with A — 180° that correspond to a non-resonant 
case become rotationally-symmetric. Note that the distributions are not radially 
uniform even in this case: the particle density is higher at the inner and outer 
edges of each ring. Mathematically, the density there gets infinitely large, because 
the radial velocity of particles vanishes in pericenters and apocenters. Figure |3] 
shows configurations of particles locked in different resonances, from 0-th to the 
3rd order. For simplicity, the eccentricity and the libration width are the same on 
all panels: e = 0.25 and A = 0.271 (36°). 



8 



Queck et al. 




-2-1012 -2-1012 -2-1012 -2-1012 -2-1012 -2-1012 

Fig. 3 Same as Fig.|2l but for objects with the same orbital eccentricity oie = 0.25 and libration 
width of A = 0.2JI (36°), locked in different resonances. From left to right: 1 : 1 (Trojans near L4), 
2:1,3:2, 4:3,5:3,7:4. 



3 Collisions 

3.1 Collision condition 

In terms of radius vectors, the collision condition is trivial: ri = r2, meaning that 
distances and true longitudes 00 + 6 of both particles coincide; 



n = r2 
©1 + 9i = ©2 



09 



(8) 
(9) 



By applying the equation of conic section and solving for (O2 and 02, we get 



(O2 = 01 -02 + (Oi 



fl2 1 - C2 



COS 02 = — 



1 

62 



fll 1 



which splits up into two solutions ("+" and "— ") 



;i+eicos0i)-l 



01 - 0, 



«2 = 

and 



CO2 = 01 - 02 



COi, 

©1, 



0+ 



arccos fcos 0^ 



02 = 2k — arccos [cos 02] . 



(10) 
(11) 

(12) 
(13) 



3.2 Relative velocity at collision 

To compute the relative velocity 1^^^^(61,(01,01 , 62,(02, 02) of two colliding par- 
ticles, we start with calculating the velocity vector of either particle. Consider a 
cartesian coordinate system centered on the star and with x-axis pointing towards 
the particle's present position, but not rotating with the radius vector. The parti- 
cle's position and velocity in this system are 

r=(;), (14) 

_fvr\_ I GM ( esin0 \ 

Denoting the velocity vectors of two colliders by Vi and V2, the fc-th power of the 
relative velocity is 

V,L = |vi-V2|^ (16) 
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Therefore, the relative velocity can be calculated by simply applying Eq. (fTSl l to 
both colliders, which makes it obvious that Vi,„p = Vimp{ei,Qi,e2,^2)- Further- 
more, at collision, 62 is determined by ei, 9i, and 62, see Eq. jlll l. Therefore, the 
impact velocity depends on three arguments only: 

Vi,„p = Vi,„p{e\,Q\,e2). (17) 

We do not give the explicit form of V,*,p, because it is rather lengthy and obtained 
by a straightforward calculation. 



4 General formalism 

4.1 Splitting of variables 

Following iKrivov et alJ (l200l l2006h . we now arrange all phase space variables 
into "useful" ones p that we keep and "dummy" ones q that we will average over. 
For our problem, we choose 

p = (e), q=((0,e). (18) 

Recall that a is a parameter, not a phase space variable. 



4.2 A-integrals 

Consider an arbitrary function defined at the collisional point, for instance the 
k-th power of the relative velocity of two colliding particles V,f„p(p],qi,P2,q2)- 
lKrivovetalJ(l2005Ll2006h have shown that the mean value of that function, aver- 
aged over all q-variables, can be expressed as 

AW(pi,p2) 

where 



A«(Pi,P2)^ 

r r (20) 
(Pi , qi , P2> q2) 5 (ri - r2) ^q(Si\)^q{ii2)diildii2- 

Hi H2 

This formula already includes the collisional condition through the factor 5(ri — r2) 
We will refer to A*^*^' as A-integral. 

With our choice of variables and setting k=\, Eq. M9\ can be rewritten as 

Vimp (e 1 , 62 ) - -7577 T , (21) 

AW(ei,e2) 

which gives the average collisional velocity between two overlapping rings of par- 
ticles: one with eccentricity e\ and another one with eccentricity 62, with angular 
elements of particles distributed in accordance with the resonance condition. 
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4.3 Meaning of A-integrals 

As was shown in iKrivov et alJ (l2005h . can be interpreted as the reciprocal of 
an "effective interaction volume". Consider again two rings formed by two subsets 
of particles with given eccentricities ei and 62- If Si and ^2 are the surface areas of 
the rings ei and 62, and 512 the area of their intersection, then the zeroth integral 
a'"' (61,62) is approximately 

The integral a'') allows direct physical interpretation, too. Since A^^' = V/mpA^"' 
^ velocity/volume, one expects that, after multiplication by the number of parti- 
cles and their collisional cross section, A^' would give the collisional rate. More 
precisely, this will be the rate of collisions between particles with eccentricity ei 
and those with eccentricity 62- This is indeed the case, see Eqs. (l28]l-(|29]l below. 



4.4 Evaluation of A-integrals 

For the actual calculation of any A-integral we insert ( fTTl l and (fTSl l into ( |20l ): 



AW(ei,e2) = 



'0)1 J(02 •'01 . 

X 5(ri -r2) 

The 5-function in ( 122b should now be expressed through orbital elements: 
8(ri-r2)=/-i5((Oi -(02)5(01 -62) , 

where the Jacobian 

/ = 



9(0)2,02) 



explicitly reads 



C2 



sin 62 1 = a' 



62(1 -4)^1 sin 921 



a 1-e^'"""" " " (l+e2Cose2)3 
Collecting all intermediate results and incorporating them into 

AW(ei,e2)= I I I I V^euQue2) 
Jm, JToy Je, Jbj 



we get 



'Ml J0>2 J<)\ 

1 (1 



62 cos 62)^ _, 



5(co2-co5)5(e2-e|) 



fl2 e2(l-e2)2|sine2| 

X (|)aj(COi,ei) (|)aj(c52,e2) iiC0iii(02 ddidd2- 



(22) 



(23) 



(24) 



(25) 



In Appendix A, this formula for the A-integral is transformed further. Three of 
the four integrals turn out to be analytically solvable, so that the final expression, 
Eq. (149b . contains only an integral over 9i, which we evaluated numerically. 
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The A-integrals, and therefore the collisional velocities and rates, have several 
useful properties, whose derivation is given in the Appendices A and B: 

1. The A-integral is symmetric with respect to its arguments: A^''\e2,e-[) = 
A«(ei,e2). 

2. The A-integral depends on the libration width A, but does not depend on the 
libration center <I>o. 

3. The limit of the A-integral at A ^ is finite, and close to the values obtained 
with small, but non-zero A. 

4. The same holds for the limit of the A-integral at ei —> 62- it is finite and is 
not far from the values calculated for close, but not equal ei and £2- 



4.5 Comparison with approach by Dell'Oro et al. 

The formalism developed h ere is based on exactly the same ideas as the one b y 
Dell'Oro and collaborators (Dell ' Pro & Paolicchil Il998l : iDeH'Oro et all Il998h . 
For instance, our A^^^-integr al ( |20] | is essentially the same as Eq. (9) or (10) in 
Ore & Paolicchil (119981) . A technical difference between the two approaches 
is that we incorporate the collisional condition directly into the integrand, through 
the function 5(r i — r2), and assume a particular functional form of the distribution 
of orbital elements (Eqs.|5]and|6ll. We then perform analytically as many integra- 
tions in the multiple integral as possible. As a result, only one integration (see Eq. 
|49] | needs to be performed numerically. The price to pay for "more analytics" in 
our approach is that it is much "heavier" mathematically, which makes a strict 3D- 
treatment impossible. Still, we deem this approach suitable for theoretical study 
of a statistical ensemble of pseudo-objects with continuous distributions of orbital 
elements. As we will see, our approach is quite convenient to explore dependence 
of the collisional velocities and rates on various parameters {e,„ax. A, p, and q). 
Besides, our approach naturally circumvents numerical difficulties that otherwise 
would arise for extreme values of these parameters (e.g., for very low or very high 
eccentricities). 

In contrast, Dell'Oro et al. use a Monte-Carlo technique to evaluate the mul- 
tiple integrals, similar to our A*^. This allows a study in three dimensions; in fact, 
they use 

p = (a,e,i) q=(n,(0,e) (26) 

(cf. Eq. [TSl l. This makes their approach ideal for study of particular collisional 
complexes in the solar system, consisting of individual objects with known orbital 
elements. Their method is particularly useful to explore effects associated with 
inclinations and longitudes of nodes. 

5 Collisional velocities 

5.1 Collisional velocity for the subsets of particles with e = ei and e = e2 

In this section, we show numerical results for the collisional velocity Vimp{e\, C2), 
calculated with the aid of (121b and (|49] |. Recall that this velocity is the average 
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V(ei,e2=0.1) |q=1, A=0.3;rl 




Fig. 4 Cuts through the contour plots of coUisional velocity Vj,„p{ei,e2 = 0.1), in units of 
^kepier, for the 2:1 resonance (solid), 3:2 (dashed), 4:3 (dotted). Top left: a shallower resonance 
with A = 0.371, bottom left: a strong one with A = O.lTt. Right panel: typical numerical results 
for V/,„p(ei ,£2 =0.1) anda2:l resonance for comparison. 

collisional velocity between two subsets of particles in the resonant family: one 
with eccentricity ei and another one with eccentricity 62- 

Figure |4] (left column) shows this velocity as a function of one of its two ar- 
guments, with the second argument fixed to 0.1. The velocity is measured in units 
of the circular Keplerian velocity, Vkepier ^ \/ GM/a. The major effect seen in 
these plots is that V/,„p(ei,e2) decreases from 62 at ei = to zero at ei = 62 and 
then increases again. The same V-shape pattern is seen in right panel, resulting 
from our test numerical integrations, in which we trapped planetesimals into a 
2: 1 -resonance with a slowly migrating, Neptune-like planet. 

The res ults were obtained using the algorith m developed by lThebault & Brahid 
(11998) and lThebauh & DoressoundiramI d2003h . 

The two left panels show the resonances of different strength, a shallow one 
at the top and a strong one at the bottom. Different curves correspond to different 
resonance numbers p and q. Although the resonant lock does affect the velocities, 
is it somewhat surprising that the effect is rather subtle. 

Fig.|5]depicts the velocity in the ei-e2-plane. The upper panels for a strong res- 
onance shows, as expected, a sharp minimum along the line of equal eccentricities 
and an extended area of low collisional velocity where both eccentricities are low, 
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Fig. 5 Contour plot of collisional velocity Vii„p(e\,e2), in units of Vitepier, for the 2:1 resonance 
(left), 3:2 (middle), 4:3 (right) with libration width A = 0. IJt (top), A = 0.3Jt (middle), and A = 7t 
(bottom, non-resonance case). Darker regions correspond to higher velocities. 



ei, 62 0. . .0.2. As far as other regions of the plane are concerned, the effects 
are highly non-linear. Particularly interesting is a sharp increase of the collisional 
velocity from moderate to large eccentricities that occurs for 3:2 and 4:3 reso- 
nances. It is due to the fact that, starting from a certain e, the resonant "clumps" 
start to overlap, see two rightmost panels in the middle row of Fig. |2] The larger 
p, the smaller the eccentricity at which the effect shows up. This is because a 
{p + q) '■ p resonance produces p clumps; the larger the number of clumps, the 
easier they start to overlap. The same effects, in a weaker form, are seen in the 
middle panels that are drawn for a shallower resonance. The velocity is highest in 
collisions between particles in highly-eccentric orbits with those in moderately- 
eccentric ones. In the non-resonant case depicted in the lowest panels, the highest 
velocity is attained in collisions between particles in highly-eccentric orbits with 
those in nearly-circular ones. The maximum possible value of V,mp(ei,C2) is \/2- 
It is achieved asymptotically when ei = and e2 — > 1- 
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5.2 Average collisional velocity in the disk 

We now calculate the average collisional velocity in the whole disk of particles. 
This is accomplished by integrating Vimp{e 1,^2) over both ei and £2 from to 
e„,ax, Eq. In such a way, we get the average collisional velocity in the disk as 
a function of the maximum possible eccentricity, Vi,np{e,nax)'- 



The additional double integration is done numerically with a Monte-Carlo method. 

Fig. m (top) shows the dependence on the libration width of the resonant argu- 
ment. The smaller A, the stronger the resonance - and its influence - gets. A = 71 
describes a non-resonant case, whereas A = corresponds to a perfect resonant 
lock. Again, the collisional velocity is affected by the resonance only weakly, but 
the effect is present. It is interesting that for moderate eccentricities, the average 
collisional velocity in a resonant belt is lower than in a similar non-resonant one, 
while for large eccentricities, the opposite is true. Similarly, middle and bottom 
panels in Fig. |6] demonstrate the dependence of Vimp on p and q. They show that 
collisional velocity does vary with the resonant integers, and may be both higher 
or lower that the non-resonant one, but the effect is usually of the order of several 
tens of percent. 

Remember that the whole treatment of the collisional velocities presented 
above has been done purely in 2D. Therefore, in real applications, the results must 
be carefully corrected for inclination terms {i Vkepier) in velocity. 

6 Collisional rates and collisional lifetimes 

Using the same methods as in the previous section, we will now investigate the 
influence of a resonance on the frequency of collisions. One might expect the 
number density of the particles in the resonant "clumps", and therefore the col- 
lisional rate, to be higher than in a similar non-resonant belt. We check whether, 
and to what extent, this expectation is true. 

6.1 Collisional rate for the subsets of particles with e = e\ and e = 62 

Like in the case of collisional velocity, we start by considering two subsets of 
particles in the resonant family: one with eccentricity ei {n\ particles) and an- 
other one with eccentricity £2 («2 particles) and "count" collisions between parti- 
cles of population 1 with those of population 2. The rate of collisions is given by 




(27) 



(iKrivov et al.Ll200ll2006l) . 



(28) 



which simplifies to 



R{ei,e2) =No^^^\ex,e2). 



(29) 
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y{emax) {p=1,q=1) 




Fig. 6 CoUisional velocity Vimpiemax) in units of V/^gp/er- Top: dependence on libration width 
A width for the 2:1 resonance. Middle: dependence on resonance parameter p for first-order 
resonances = l,A = 0.l3i. Bottom: dependence on the order of resonance^ for p = 3, A = 0.l7i. 
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Here, a is the collisional cross-section for equal-sized particles, and the front fac- 
tor A'^ depends on what we mean by "collisional rate". If R{ei,e2) is the number 
of collisions per unit time that a particle of population 1 has with any particle of 
population 2, then = «2- Conversely, if we consider a particle of population 2 
colliding with population 1 particles, then N = n\. Finally, to get the total number 
of collisions between particles of both families occurring per unit time, we should 
setN = ni«2. 

However, the above formulas are two-dimensional. For instance, a should be 
understood as 2s, where i is the radius of equal-sized particles. Any physically 
meaningful calculation of collisional frequencies requires a 3D treatment. We thus 
introduce an approximation for a 3D A-integral, 

A»(.)^^, (30) 

where /j = 2r sine w la^es sine is the disk thickness at the annulus location, e being 
the half-opening angle of the disk. Thus the 3D collision rate is 

R(eue2)= No A'ijlie, €2), (31) 



9 (1 

now with the "usual" collisional cross section a = Ins . We note that A,, 



3D 



I.e. 



R without the fa ctor No, is exactly what is usually called "intrinsic collisional 
probability" (e.g. lGreenberg.ll982l : Davis & Farinellaill997t) . 

Instead of dealing with the full collisional rate ( [3TT l. it is convenient to intro- 
duce a dimensionless factor ^(ei, 62): 

R{ei,e2)=R{ei,e2)/Ro, (32) 

where 

j;0= ^ ^^^kepler ^^^^ 

2n^ h al„. 



res 



is the approximate "particle-in-a-box" value of the collisional rat^E In the first 
order, /?o does not depend on eccentricities of the colliding particles and depends 
on their inclinations through the divisor h only. It describes the rate of collisions 
in a non-resonant, and therefore rotationally-symmetric, ring of objects. Again, 
Rq/{No) is the intrinsic collisional probability in such a ring. Therefore the devi- 
ations of the dimensionless collisional rate R from unity would describe effects of 
the resonance, as well as some corrections due to non-zero eccentricities. 

Figure I2] (left column) shows the dimensionless collisional rate ^(ei,e2) as 
a function of one of its two arguments, with the second argument fixed to 0.1. 
For large A, i.e. for a weak resonance, shown in the top left panel, the collisional 
rate is almost independent of ei, being close to the "particle-in-a-box" value. For 
stronger resonances (bottom left panel) R{e\,e2) peaks at intermediate values of 
eccentricity ei. The stronger the resonance, the more pronounced the maximum. 
For the 4:3 resonance and A = O.lTi, the maximum collisional rate at ei ^ 0.2 is 
about 2.5 times larger than the non-resonant rate. Interestingly, the larger the reso- 
nant integer p, the smaller the value of eccentricity, at which the collisional rate is 

' See, e.g., Eq. (18) in lKrivov et aLH2007h . Note that in Eqs. (15), (17), and (18) of that paper, 
an additional factor 1/2 is erroneously missing. 
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Fig. 7 Cuts through the contour plots of collisional rate R{ei,e2 = 0.1), for the 2:1 resonance 
(solid), 3:2 (dashed), 4:3 (dotted). Top left: a shallower resonance with A = 0.3%, bottom left: 
a strong one with A = O.Iti. Right panel: typical numerical results for R{e\,e2 = 0.1) in a 2:1 
resonance for comparison. 



the highest. Another finding is that, for sufficiently large eccentricities, resonances 
may decrease the frequency of collisions. For comparison, in the right panel we 
show a typical result of our numerical integrations, in which we trapped planetes- 
imals into a 2:1 -resonance with a slowly migrating planet. Notwithstanding large 
oscillations, caused by a limited number of particles integrated, it shows the same 
tendency of the collisional rate to slightly increase with e\. 



Fig. [8] depicts the collisional rate in the ei-e2-plane. The upper and middle 
panels show that collisions are most frequent between particles, whose orbital 
eccentricities are moderate and not very different from each other. Interestingly, 
these are exactly the eccentricities for which collisional rates in the non-resonant 
case are the lowest (bottom panels). For some resonances, another region of higher 
collisional rates is observed at very high eccentricities. This effect has the same 
origin as a similar effect in the collisional velocities. Fig. [5] 
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Fig. 8 Contour plot of the dimensionless coUisional rate R{e\ , et), for the 2: 1 resonance (left), 
3:2 (middle), 4:3 (right) with libration width A = O.lTt (top), A = 0.3n (middle), and A = n 
(bottom, non-resonance case). The darker the contours, the higher the rate. 

6.2 Average collisional rate in the disk 

Like in the case of collisional velocity, we now average over both eccentricities in 
order to obtain the collisional rate in the whole disk containing n objects: 

)iSfe{ei;emax)de2dei. 

J ei J €2 

Again, the front factor A'^ depends on what we intend to describe with the function 
Ri^max)- If Ri^max) is the number of collisions per unit time that a certain particle 
may have with all other particles in the belt, then N = n.To get the total number 
of all collisions occurring in the ring per unit time, we should set N = n^/l. 

Figure |9] shows the numerical results. The top panel focuses on the influence 
of the libration width A, whereas the middle and bottom ones illustrate the depen- 
dence on p and q. It is seen that R{e,„ax) has qualitatively the same properties as 
R{e\,e2)'. independence of e,„ax out of resonance and a maximum at intermediate 




Fig. 9 Dimensionless coUisional rate R(e,„ax) in the disk, Top: dependence on libration width 
A width for the 2:1 resonance. Middle: dependence on resonance parameter p for first-order 
resonances l,A = 0.l3i. Bottom: dependence on the order of resonance^ for p = 3, A = 0.l7i. 
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eccentricities for deeper resonances. The maxima shift to shghtly lower eccen- 
tricities, when either the libration width decreases, p increases, or the order of 
resonance q decreases. All the curves are generally flatter than the non-averaged 
ones. 

The properties described above can also be understood analytically. Consider, 
for instance, the eccentricity which yields the highest collisional rate, i.e. the 
position of the maxima in Fig. |9] By calculating the limit of the A^'' -integral at 
A ^ and ei,e2^ e (see Appendix), one gets an approximate equation 



p + q 



(l-e,)3/2 



(35) 



For example, {p,q) = (1,1) gives = 0.31; {p,q) = (2,1) gives e^, — 0.19; 
{p,q) = (3,1) gives e^, = 0.14; {p,q) = (3,2) gives e^, = 0.24. A comparison with 
the positions of relevant maxima in Fig.|5l calculated for a finite A, shows that the 
former lie slightly to the left from the former, as expected. 

The collisional rates are intimately connected to collisional (life)times of the 
disk particles. The average collisional lifetime of an object in the disk is 



)^e{euemax)de2dei 



Li Ie2 ^e{e2;e,„ax)i^e{e\\emca)de2dei 



(36) 



or simply 



T^coui^^max) 



(37) 



where R{e,„ax) is given by Eq. 041) with N = n. 

Remember that we have calculated R in quasi-3D. More exactly: in the colli- 
sional rates, which are velocity / volume, we do take into account the inclina- 
tions when calculating the volume, but still ignore those in velocity. 

Non-inclusion of the term ^ / Vkepier in velocity introduces an error of the same 
order of magnitude as in the case of collisiona l velocities consider ed in Section 5 
(see discussion around Eq. (A14) and (A15) in lKrivov et aU[2005) . 



7 The 1:1 resonance: Trojans 

We now discuss a special case ^ = 0, or 1:1 resonance, which corresponds to a 
Trojan cloud of objects at one of the trigonal Lagrangian points. All the formulas 
derived and discussed above are valid in this case, but the results reveal impor- 
tant qualitative differences from the first- and higher-order resonances. Figure [TOl 
presents the collisional velocity and the collisional rates for Trojan clouds with 
different maximum eccentricities e^ax and different libration amplitudes A. 

Like in the case ^ > 0, V,„,p vanishes when e,„„^ goes to zero. However, unlike 
for other resonances, R{e,nax) has a maximum at e,„ax = 0. What is more, that 
maximum collisional rate goes to infinity when A ^ 0. Mathematically, this is 
explained in Appendix |b1 for instance, Eq. ( l59b is no longer valid. Geometrically, 
the explanation is obvious (see Fig.[3]|. For q > I, even in the case where e = 
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Fig. 10 The case of Trojans. Top: collisional velocity Vimp{^max) in units ofVkepiey^ cf. Fig.Hl 
Bottom: dimensionless collisional rate R(emax)', cf. Fig.|9] Different curves correspond to differ- 
ent libration widths A. 
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and A = the objects form an (infinitesimally narrow) ring. For ^ = 0, in that 
case the cloud simply shrinks to the Lagrangian point — in other words, all the 
objects reside exactly at one and the same point and have zero relative velocities. 
Their volume density is infinitesimally high, and so is the collisional rate. For 
small non-zero f^max and A, frequent collisions are expected. In contrast to this, 
the velocity would be lower, but the effect here is almost negligible, because the 
velocities would be dominated by 3D-terms coming from inclinations, which we 
do not consider in our model. The implications for the Trojan asteroids in the solar 
system are discussed in the next Section. 
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8 Applications 

The results may have a variety of astrophysical implications, whose detailed anal- 
ysis is beyond the scope of this paper. Here we consider some of them and only 
briefly. 

1. Resonant families of small bodies in the solar system. Our solar system is 
known to possess a number of small body populations locked in mean-motion 
resonances. These are Trojan swarms of Jupiter and Trojans of other planets, var- 
ious asteroidal families in the main belt, Plutinos and Twotinos in the Edgeworth- 
Kuiper belt. We note again that the "extreme" analytics that we employed here 
required a number of simplifying assumptions. For particular resonant families 
in the solar system, our model is obviously oversimplified to provide quanti- 
tatively accurate results. Indeed, each of the families counts typically tens to 
hundreds of known objects, whose orbital elements are known at the individ- 
ual level. Therefore, to study their co l lision al evolution, direct A'^-body integra- 
tion (e.g. iThebault & Doressoundiraml l2003h or ano ther version of statistical ap- 
proach, developed by Dell'Oro and collaborators re.g. lDell'Oro & Paolicchilfl998l) . 
should be preferred. Nevertheless, we make some qui ck comparisons of our results 
with those obtained with more accurate methods (Ma rzari et al.lll99^ : lDeirOro et al.l 
Il998l) . In principle, we could choose any resonant group, as long as the number of 
known members is sufficiently high to justify application of statistical approach. 
We decide to consider Trojans, because for them, as noted above, the effects of 
resonance on the collisional velocities and rates are the strongest. 

Assuming for either of the two Jupiter Trojan swarms e„^ax = 0.1 and A — 
O.Iti w 20°, Fig. 1101 suggests R w 10, implying that the collisions are an order 
of magnitude more frequent than simple particle-in-a-box estimates would pre- 
dict. This result can be compared to the collisional rates among (non-resonant) 
main-belt asteroids. If the latter had the same semimajor axes and the same distri- 
bution of inclinations as Trojans, the rates would be 10 times lower. According to 
Eq. ( 1331) . however , the collisional rate scales with the semimajor axis and mean 
inclination (which enters the rate through h or e in Eq.l33b as 7? oc a~''^^ {i)~^ . The 
main-belt asteroids are roughly twice closer to the Sun than the Trojans, and their 
mean inclinations of (/) = 8° are twice smaller than that of Trojans (15°). This 
would give a factor of 20 lower rates, which would nearly compensate the ten 
times higher R. Therefore, the expected collisional rate for Trojans, despite their 
larger heliocentric distance and broader vertical distribution, should be compara- 
ble with that for asteroids in the main belt. Note that the result is rather sensitive 
to the choice of e„,ax and A and that a more realistic distribution of eccentricities 
and resonant arguments can change the results signific antly. With this caveat, our 
estimate is in a reasonable agreement with a result by iMarzari et alJ (fl996.') who 
found that the intrinsic collisional probabilities for Trojans are about twice as high 
as in the main asteroidal belt. Unlike the collision rate, the effect of the resonance 
on the collisional velocity is minor. In fact, given rather large inclinations of Tro- 
jans, Vimp will be dominated by inclination terms on the order iVkepier that are 
ignored in our treatment. As a result, we expect nearly the same impact velocities 
for Trojans as for the main-belt asteroids, and, as a co nsequence, the same c olli- 
sional outcomes. Th is fully agrees with conclusions of iMarzari et alJ (Il996h and 
iDeirOroetai] (fl99kl as well. 
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Similar estimates can, of course, be made for the collisional evolution of Kuiper 
belt families. For eccentricities lower than e < 0.3, and the amplitude of libration 
of the resonant argument of tens of degrees, the collisional rates may be up to a fac- 
tor of two higher than for non-resonant populations. The collisional velocities will 
not be appreciably affected by the resonarices. F i nally, applications to asteroidal 
families are also possible ^cf . iMarzari et aliri995l : lDeirOro et al l l200lLl2002h .In 
this case, however, the theory developed in this paper needs to be adapted to the 
case of internal resonances. 

2. Clumpy debris disks. Our model may be particularly useful in the cases 
where orbital distributions in the small body families are poorly known or/and 
where the number of objects in the family is far too large for A'^-body integrations 
to work. The first fully applies to extrasolar "asteroid belts" around other stars, 
whose existence is suggested by indirect, and very uncertain, methods. The sec- 
ond applies directly to dust clouds such belts should produce. This makes debris 
disks an ideal application for the model, as it enables easy estimates for different 
central stars, parameters of presumed planets and pl anetesimal be l ts, and various 
resonances in which planetesimals might be locked. iKrivov et al.l (120071") consid- 
ered two possible scenarios for formation of clumps observed in resolved debris 
disks. In a standard scenario, the Poynting-Robertson force delivers dust from 
outer regions of the disk to locations of external mean-motion planetary reso- 
nances w ith a planet; dust gra ins get captured and form visible clumps. In another 
scenario (I Wyatd I2003L 120061) . a population of invisible planetesimals resides in a 
resonance with the planet, such as Plutinos and Trojans in the solar system; the 
dust produced by these bodies would stay locked in the same resonance, creating 
the dusty clumps. With the help of simple analytic models, iKrivov et al.l (l2007l) 
showed that the first scenario works for disks with the pole-on optical depths be- 
low about ^ 10^^-10"^. Above this optical depth level, the first scenario will 
generate a narrow resonant ring with a hardly visible azimuthal structure, rather 
than clumps. The model of the first scenario was based on the particle-in-a-box 
formulas for collisional rates, ^ = 1 . The fact that the corrections to R found here 
for the resonant case are only moderate lends further support to the conclusions 
of that paper. On the other hand, these corrections may lead to some quantitative 
changes. The "critical" optical depth, or fractional luminosity, of a debris disk in 
which the standard scenario may be efficient, shifts towards values that are by a 
factor of several lower. 

9 Conclusions 

The subject of this paper is "statistical" celestial mechanics of the restricted three- 
body problem. Specifically, we have investigated the averaged collisional veloci- 
ties and rates of collisions in an ensemble of planetesimals orbiting a star, locked 
ina{p + q) : p-resonance with a planetary perturber that revolves around the same 
primary, and having a distribution of eccentricities. 

We use a statistical approach, ba sed on the same philosophy as the one devel- 
oped by Dell'Oro and collaborators (iDell 'Oro & Paolic"ctii[ll998tlDeir Oro et al.L 
[l998). The main difference between the two approaches is that we treat a statisti- 
cal ensemble of pseudo-objects with continuous distributions of orbital elements, 
rather than a finite set of objects with given orbital elements. 
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Our main findings are as follows: 

(i) The mean collisional velocity is nearly the same as in the case when the 
family is non-resonant. In other words, it does not differ notably from the one pre- 
dicted by the particle-in-a-box method: i.e., that velocity increases almost linearly 
with the maximum eccentricity e,„ax of objects in the resonant family. 

(ii) If the eccentricities of the objects are low to moderate, the mean colli- 
sional rates are higher within a resonance than out of it. The enhancement in the 
collisional frequency increases from shallow resonances to deeper ones, i.e. with 
decreasing libration amplitude of the resonant argument. For the libration widths 
of tens of degrees, the collisional rate may increase by about a factor of two. There 
is a nonlinear dependence of the collisional rate on e,„aA : the rate is highest at a 
certain e„,ax- The smaller the resonant integer p or the lower the order of resonance 
q, the smaller that emax- The maximum collisional rate is highest for the first-order 
resonances, q = I. 

(iii) In the special case of the primary, 1 : 1 resonance, the effect is the largest. 
The collisional rates increase, and the velocities decrease, with decreasing eccen- 
tricities and/or libration amplitudes in the family. 

(iv) All the conclusions presented above apply to collisional velocities and 
rates averaged over the swarm of objects. Therefore, the statement that rates and 
especially velocities of collisions are only moderately affected by the resonant 
lock should not be misinterpreted. Studies with A^-body codes show that locally, 
namely at longitudes that correspond to the location of the clumps, both velocities 
and rates may be orde rs of magnitude higher that outside the clumps (see, e.g.. 
Fig. 4 in Wyat]| |2006|) . Such studies show, furthermore, that clumps are indeed 
regions where most of the collisions take place and where most of the collisional 
debris emanate from, etc. In contrast, our results are more pertinent to the global 
and long-term evolution of the whole resonant population, e.g. collisional lifetime 
of objects, mass outflow from the system and so on. To illustrate the difference, 
let us discuss the lifetime of individual object. It is not possible to say whether an 
individual planetesimal "belongs to a clump" or "it does not". It spends part of its 
time in the clump and part out of it. Averaging over the orbits is needed to predict 
its fate statistically, and that is exactly what we have done in this paper. 

The results described above rest upon a number of simplifying assumptions 
that were made to keep the problem analytically tractable. An incomplete list of 
them includes circular planetary orbit, treatment in the 2D, like-sized objects, tak- 
ing the same semimajor axis for all orbits, uniform and constant distribution of 
eccentricities and resonant arguments, a single isolated resonance without interac- 
tion with background non-resonant objects and other resonant families. However, 
even with all those assumptions, the mathematical complexity of the derivations 
is rather high, so that the prospects to relaxing most of them are questionable. 
From this viewpoint, we have demonstrated that an analytic approach has severe 
limitations, both regard to the accuracy of results and their applicability range. 
Nevertheless, our model has a useful genericity, allowing one to quickly estimate 
the strength of the resonant effects on the collisional evolution for various reso- 
nances, different libration amplitudes, different dispersions of eccentricities, etc. 
Last but not least, as with any analytic study of a simplified problem, it offers a 
clear dynamical explanation of why and how the resonant lock alters the rates and 
velocities of collisions between the objects. 
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A Transformation of A-integral 

Here, we will transform Eq. i25\ to the form suitable for efficient numerical calculations. 

Two out of four integrations in Eq. i25i . those over Wi and 82, can be done immediately with 
the help of the two-branch collision condition. Since cos 9 J = cos 62 = cos 9^ and | sin 9^ | = 
btain: 



I sin 92 I = I sin 92 1 we obtain: 



1 /•2''^g^(l+e2COS0|)3 



a2e2il-e2^y Jo |sin95| 

X llvL i'^uQuei) / da5,(^cj(coi,9i)(|)oj(col,9|). (38) 

^ Jo 

From now on, the superscript "c" of and 01 will be omitted. We stress, however, that both 
variables are functions of ej, coi and 9|, and 62, as calculated from Eqs. i fTot and <1 lb . The 
A-integral takes the form 



A«(^i,e2) 

(l_e2)l/2(l_, 2)1/2 .2. 



d9i 



X 

± 



167t2A2a2 Jo |e2sin92|(l-l-eicos9|) 

, rliz 

(ei,0,,6.2) / dM, //1//2 (39) 



where the Heaviside functions 

Hj=H[^o-A<<i>j <<t>o+A], j = (l,2) (40) 

describe libration of resonant arguments 'I>i and <t>2 around <I>o (e.g. 71 or 0). 

The integral / iS5\H\H2 is analytically solvable. Resonant arguments and $2 can be 
written as 

Hj=H[-A<{pmi+Xj}2K<A\, i=(l,2), (41) 

where {■}2n denotes a modulo 2ji operation which returns values between [— Jt, Jt]. E.g. {0}2n = 
0, {7t/2}2n = It/2, {— 7t/2}2n = {37i/2}2k = — 7t/2, etc. The arguments in the above equations 
are 

%\ = {p + q)M(ei,Q\) - <fo, (42) 
X2 = (p + q)M{e2,Q2) - tfo - p(02-9i). (43) 

Here, M = M{e, 9) is the mean anomaly given by Eqs. ||3}-||4)- 

Let us first consider the integral Jq^ dWiHi . Here, the offset X\ is unimportant and it is easy 
to see that the integral evaluates to 

t fin 

dMi//i = / dWiHl-A < {pmi}2n <A] = 2A. (44) 

Jo 

Similarly, / dc5i//|//2 can be simplified to 

/ dmiHiHT = / dmi H[-A<{pWl}^,t<A] 
Jo Jo 

X H[-A<{pTa^+Ax}2^<A], (45) 

while the shift Ax = {X2 - Zi }2k 

AX = {{p + q) [M(e2,92)-M(ei,0i)]-p(92-9i)}2„. (46) 
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This can be further simplified by noting that we are "overlapping" intervals [—A, A] and [A^ — 
A,Ax+A] in modulo 27i space. Therefore 



drai//|//2 = 2A 1 



|Ax| ^ Jt 
2A - A' 



2A 



H 



2A 



(47) 



For A < Ji/2 this branch is always taken. Another branch is 



dMi//i//2 = 2A (2 - -V 

\ AJ 



|Axl ^ 7t 
2A - A' 



(48) 



Note that |A%| < Jt and A < Jt. In the limiting case A = Jt we obtain / A(ii\HiH2 = 2A. 
The final result for the first branch (always taken for A < Ji/2) is: 



A«0 



(l_e2)l/2(i_g2)l/2 



87t2Aa2 



Jo 



de 



1 



62 sin 02 1(1 +ei cosBi) 
k, ^ J, |Ax±| 



^V^n;{euQue2)[l 



2A 



X H 



|AX^ 
2A 



(49) 



and similar to the second branch. The only remaining integral is evaluated numerically using a 
Monte-Carlo method. 

The above calculation assumes that the distribution of the resonant argument within the 
libration band is uniform, Eq. (|5j- If necessary, they can be repeated with a more realistic 
distribution. It is sufficient to replace the Heaviside function in Eq. Js) by another function 
/({O — <I>o}2jt)i for instance a trigonometric one or a Gaussian, and re-do the calculation de- 
scribed here. 

A more detailed analysis of the above equations allows one to establish several useful prop- 
erties of the A-integrals. 

First, it is symmetric with respect to its arguments. Although it is not quite evident from 
l |49l l, one can use Eq. (TT) to change the integration variable from 9) to 82. The result will be 
identical to A(e2.ei), and therefore K(e2,e\) = A(ei,e2)- Of course, we have also tested the 
symmetry numerically. 

Second, <49l l readily shows that the simplified A integral depends on the libration width A, 
but does not depend on the libration center <l>o. 

Third, it is possible to calculate limiting cases of Eq. <49t and to get an idea of how large 
the effect of resonance on collisional velocities and rates could be (Appendix [bJ. 



B Limiting cases of A-integral 

Here we calculate the double limit of the a'*~'' {e\ , e2)-integral at A — > and at e\ — > ^2- The first 
corresponds to an "exact" resonance. The second "emulates" averaging over both eccentricities 
in A(ei,e2)- 

We start with the limit A ^ of Eq. l l47t : 

/•271 

lim / dcoi//i//2 =4A2 5(Ax), (50) 
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which completely eliminates the A dependence, as Eq. J39t contains A^^. Thus the A—>0 limit of 
the A-integral is finite. Furthermore, the 5-function immediately resolves the remaining integral 
over 01 to an explicit expression since 



5(Ax(ei))=^ 



aAx(ei 



S(ei-ef'), (51) 

where 9','' are the roots of the equation Ax(6*i'') = 0. In practice, the calculation of Eq. <51t is 
cumbersome, as it requires inversion of the Kepler equation. 

We now consider the ei ^ ei limit, confining our analysis to the case A(*='). After some 
algebraic manipulation, it can be shown that 

ei^ei |e2SmG2 I 

lim ^""P' 
|e2sin92 

Here, the "±" branches of the collision condition have slightly different meaning than in (fT2t- 
i fTst . With the superscript "+" we denote the case 92^^ 9i, and with "— " we denote 92*^ — ► — 9i. 

For the remaining nonzero "— " branch, the condition Ax(9j'') = is satisfied at 9j'' = and 
9p' = 71. The integral A''' is again finite, except at points f i = £2 = when 



lim I ^""''^ _| = Vkepler- (53) 



3A3c(9i 



39] 



= 0. (54) 



At points ej = 62 = Eq- HD implies that A''' is infinite, and those values correspond to 
maxima in Fig. [8] Note that the maximum rates, R(e,t,e*), will rise to infinity when A —> 0. 
However, the A ^ limit of A*'' (ei , ej) integrated over ei and ej, will stay finite. As a result, 
the maxima in Fig. [9] will not be infinitely high when A—>0. 
Equation <54l l reduces to 



301 



(55) 



The case with 0^^' = n does not have a solution, while for 0,'' = the equation takes the form 



{i-e,)y^ p ■ 

The equation Ax{Q^i ') = usually has more than two roots (e.g. for ei = ^2 = those are 

0[''' = ji^/q, where j is integer and \j\ < q). In certain cases these additional roots can also 
simultaneously satisfy Eq. \5A\ leading to additional maxima (see Fig.[8]middle and right). 

In order to find other maxima of the A-integral, it is easier to go back to Eq. {55). For q>\ 
the only solution to Eq. ( |55t is a pair (0i = Q\d,e\ = £2 = ^rf). where 

^ ^ {^-e/?l'^(p + q)lp + e,-\ ^ ^^^^ 
led 

Q\d = ±2arccos [v^] . (58) 

It can be proved that x > and dx/ded < (for q/p > and < ed < I), which has a conse- 
quence that ed = (when x = I and 0i^ = 0) is the smallest ed for which the A-integral has a 
maximum. Then, the other maxima (if any) can be obtained using Eq. iSSi to check which pairs 
(Qirfi^i/) satisfy A^ = 0. The density plot of |Ax(9i,e2 = ei)|, with the lines (9i = 6id,ei = 



Collisions in resonant belts 



29 




Fig. 11 The 2D plot of |Ax(ei,e2 = From left to right: 2:1, 3:2, and 4:3 resonance. Black 
stripes in the density plot correspond to Ax = 0. The intersection between the black stripes and 
the solid line (ed,6irf) are the maxima. The white dot marks the first maximum, (e*,0). 



^2 = ^d) overplotted, is given in Fig.[TT]for the same three resonances as in Figs.|5]and[8] 2:1, 
3:2 and 4:3. 

We finally note that for g > 1 , 

lim a(i) = (59) 

All of the equations in these Appendices hold for the case ij' = (1:1 resonance, Trojans) as 
well, except when an equation has a explicit 1/g. In particular, g = in Eq. \56\ leads to e^, = 0, 
which agrees with Fig.flOl 

So far, we considered the integral A(*=''. Of course, the integral A'*^''' is not less important, 
because it is needed to calculate the coUisional velocity, see. Eq. M\\ . For k = the "+" branch 
(Eq. I52t diverges when ei — > ^2- Therefore, a'*^''' also diverges, and the coUisional velocity 
l l2l1 i goes to zero, as it should (see Fig.|4}. 
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